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Abstract. Let N{f) be a number of nodal domains of a random 
Gaussian spherical harmonic / of degree n. We prove that as n 
grows to infinity, the mean of N{f)/n'^ tends to a positive constant 
a, and that N{f)/n^ exponentially concentrates around a. 

This result is consistent with predictions made by Bogomolny 
and Schmit using a percolation-like model for nodal domains of 
random Gaussian plane waves. 



1. Introduction 

Let Tin be the 2n + 1- dimensional real Hilbert space of spherical 
harmonics of degree n on the 2-dimensional unit sphere equipped 
with the norm. For / G Hn, put Z{f) = {x G : f{x) = 0}. 

Let A^(/) be the number of connected components of Z{f). The famous 
Courant nodal domain theorem [31 Chapter VI, § 6] states that A^(/) ^ 
(n+ 1)^ for all / G Hn- On the other hand, H. Lewy [5] showed that no 
non-trivial lower bound is possible: one can find spherical harmonics 
/ of arbitrarily large degree with A^(/) ^ 3. The question we want to 
discuss here is: What is the "typical" value of N{f) when the degree n 
is largel To give the word "typical" a precise meaning, let us consider 
the random spherical harmonic 



f=Y. ^^^^ 



where are independent identically distributed Gaussian random vari- 
ables with E^^ = and {Yk} is an orthonormal basis of Tin, so 
^11/11 r2/c2N — 1- is hard to see that / (as a random function) 

L (s ) 

does not depend on the choice of the basis {Y^} in Tin- 

The same question can be raised in other instances of smooth ran- 
dom functions of several real variables, e.g., for random trigonometric 
polynomials of large degree n. We are not aware of any rigorous treat- 
ment of this question, though we know two encouraging attempts to 
tackle it in very different contexts. In the paper [7] (motivated by 

1 



2 



FEDOR NAZAROV AND MIKHAIL SODIN 



some engineering problems), Swerling estimated from below and from 
above the mean number of connected components of the level lines 
Z{t, f) = {f = t} of a random Gaussian trigonometric polynomial / 
of two variables of given degree n. His method is based on estimates 
of the integral curvature of the level line Z{t,f). The estimates are 
rather good when the level t is separated from zero, but as t — they 
are getting worse and, unfortunately, give nothing when t = 0. 

A few years ago Blum, Gnutzmann, and Smilansky [Tj raised a ques- 
tion about the distribution of the number of nodal domains of high- 
energy eigenfunctions. In the ergodic case, in accordance with Berry's 
"random wave conjecture" , they suggested to find this distribution for 
Gaussian random plane waves and performed the corresponding numer- 
ics. To compute this distribution, Bogomolny and Schmit suggested in 
an elegant percolation-like lattice model for description of nodal do- 
mains of random Gaussian plane waves. It agrees well with numerics, 
but completely ignores the correlation between values of the random 
function / at different points, and apparently it will be very difficult 
to make it rigorous. 

In this note, we will show that, in accordance with one of the Bo- 
gomolny and Schmit predictions, E,N{f)/n'^ tends to a positive limit a 
when n oo. Moreover, we show that the random variable N{f)/n'^ 
exponentially concentrates around a: 

Theorem 1.1. There exists a constant a > such that, for every 
e > 0, we have 

where c{e) and C{e) are some positive constants depending on e only. 

Remark 1.2 (Sharpness of Theorem ll.il) . The exponential decay in n 
in Theorem 11.11 cannot be improved: in Section [H] we show that, given 
a positive and arbitrarily small n, P|A^(/) < kv?' } ^ e-^^")". On the 
other hand, our proof of Theorem 1 1 . 1 1 gives a very small value c{e) > e^^ 
and it would be nice to reduce the power 15 of e to something more 
reasonable. 

Remark 1.3. The model proposed by Bogomolny and Schmit also 
predicts that the variance of the random variable A^(/) grows with n 
as hn^ with some constant 6 > 0. 

Remark 1.4. For any spherical harmonic f & Ti, the total length of 
its nodal set Z{f) does not exceed Const n. Therefore, Theorem 11.11 
yields that, for a typical spherical harmonic, most of its nodal domains 
have diameters comparable to 1/n. 
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The proof of Theorem 11.11 goes as follows: 

I. First, we prove the lower bound EA^(/) ^ const n^. This part of the 
proof is rather straightforward and short. 

II. Then we prove the exponential concentration of the random variable 
N{f)/n^ around its median. This part is based on two ingredients: 

(i) the uniform lower continuity of the functional / ^— A^(/) with 
respect to the L^-norm outside of an exceptional set E G H of expo- 
nentially small measure; 

(ii) Levy's concentration of measure principle. 

III. In the third part, we prove existence of the limit lim EA^(/)/n^. 

n— >oo 

In this part, we use existence of the scaling limit for the covariance 
function E|/(a;)/(y)}. 

Note that in the proof of Theorem 11.11 we use only relatively sim- 
ple tools from the classical analysis, which we believe may work in 
a more general setting of random functions of several real variables, 
while it seems that the Bogomolny-Schmit model is essentially a two- 
dimensional one. 

Notation. Throughout the paper, we denote by c and C positive nu- 
merical constants whose values may vary from line to line. As usual, 
the constants denoted by C are large, while the ones denoted by c are 
small. In the cases when we need to fix the value of some constant, we 
assign to it a certain index, for instance, Cq and Cq. 

By T>{x, r) we denote the spherical disk of radius r centered at x, by 
D{y,t) we denote the Euclidean disk of radius t centered at y. 

By a we denote the spherical area measure with normalization cr(S^) = 
1, and by m we denote the (Euclidean) area measure on the plane. 

By II ■ II we always mean the L^(§^)-norm. 

Given a set K, we denote by K^^ the ^-neighbourhood of K. We 
apply this notation both to subsets of 7i„ and the L^-distance, and to 
subsets of and the usual spherical distance. 

Notation A < B and A> B means that there exist positive numeri- 
cal constants C and c such that A ^ C ■ B and A ^ c- B. If A < B and 
A> B simultaneously, then we write A ~ 5. Notation A <^ B stands 
for "much less" and means that A ^ c - B with a very small positive c; 
similarly, A ^ B stands for "much larger" and means that A ^ C ■ B 
with a very large positive C. 

Acknowledgments. We learned about the problem considered in this 
note and about the works [H [2] from Zeev Rudnick. We thank him as 
well as Leonid Polterovich, Boris Tsirelson and Steve Zelditch for very 
helpful discussions. 
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2. Main tools 

2.1. Spherical harmonics. We shall need a few standard facts about 
spherical harmonics of degree n. Most of them can be derived either 
from the fact that they are eigenfunctions of the Laplacian on the 

sphere corresponding to the eigenvalue n{n + 1) or from the fact that 
they are traces of homogeneous harmonic polynomials of degree n on 
the unit sphere. Everywhere below we assume that n ^ 1. 

Claim 2.1. The scalar product in the Hilbert space Tin is invariant 
under rotations (and reflections) of the unit sphere. As a consequence, 
the distribution of the random spherical harmonic f is also rotation 
invariant. 

Claim 2.2. For any f G Hn and any point x e & , we have 

\m\'<n'f f; 

Jv{x,l/n) 

|V/(x)r<n^/ f; 

JV{x,\/n) 

|VV/(a;)r<n6/ f\ 

Jv{x,l/n) 

Claim 2.3 (Length estimate). For any f e 7i„ that is not identically 
0, the total length of Z{f) does not exceed Cn. 

The next claim follows from the classical Faber-Krahn inequality: 

Claim 2.4 (Area estimate). For any connected component Vt of E!^ \ 
Z(f), we have Area(Q) > n~^. 

Next, we bring several classical facts about the Legendre polynomials 

Pjx) = -^^^ ^. Note that PJl) = 1. 

Claim 2.5. The function Yo(ip,9) — \j2n + lP„(cos^) is a spherical 
harmonic of degree n with \\Yq\\ — 1. Here {(p,9) are the longitude and 
the co-latitude on the sphere 

The function Yq is called a zonal spherical harmonic. 

Claim 2.6. There exists a positive constant Cq such that 

everywhere on 
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This claim is a combination of two classical results: 

(i) The function x ^— n{n + 1)P^ + (1 — x^)(P^)^ increases on [0, 1] 
(this is Sonine's theorem [8| Section 7.3]). 

(ii) |P2m(0)| = Qm and |P2m-i(0)| = 2m5(m, where = ■ 
This follows, for instance, from the recurrence relations 

(n + l)P„+i(x) = {2n + l)xPn{x) - nP„_i(x) 

and 

FU,ix) = xP:Xx) + in + l)P^ix) 
[H Chapter VII, § 3]. Recall that by Wallis' formula, gm ^ 

Claim 2.7 (Reproducing kernel in Tin)- 



V Y,{x)Y,{y) = P„(cose(x,y)) , 

k=—n 

where Q{x,y) is the angle between the vectors x,y & S^. In particular, 

n 

Y,'{x) = 2n+1. 



k=—n 



The next two facts can be found in Szego's book [8j (Theorems 6.21.2 
and 8.21.6 correspondingly). 

Claim 2.8. Suppose x^ = cosO^ are zeroes of Pn enumerated in decay- 
ing order: 

+1 > xi > ... > Xn> -1, < 6'i < ... < 6n < re. 

Then 

2z/ - 1 ^ 2u 

71 < Oy < TT, = 1, 2, ... n . 

2n + l 2n + l ' 

Claim 2.9 (Hilb's asymptotics). 

1 /2 

P„(cos^^)= (-^\ J^{{n + \)e)+R, 
\smt/ J 2 

where 

' e^/^0{n-^/^), C/n^9<: n/2 
and Jo is the zeroth Bessel function. 



R 



Note, that we shall use Claim [2^8] only for v = 1 and 2, and Claim [2^ 
for ^ ^ ^ C/n. 
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2.2. Probabilistic claims. We shall also need a few classical facts 
about the Gaussian random vectors in spaces of high dimension. 

Claim 2.10 (Bernstein's concentration of norm). 

P{||/|| >2}^e-^". 

The next result follows from the Gaussian isoperimetric lemma which 
is due to Sudakov-Tsirelson [6J and Borell [3]: 

Claim 2.11 (Levy's concentration of Gaussian measure). Let F C Tin 

be any measurable set of spherical harmonics. Suppose that the set F^p 
satisfies F{F+p) < f . Then P(F) ^ 26"'=^'". 

Claim 2.12 (Independence of / and V/). If x & then f{x) and 
V/(x) are independent Gaussian random variables. Also, due to rota- 
tion invariance, we can say that E|/(x)p = 1, E|V/(x)p < , and 
that the distribution of Vf{x) is rotation invariant on the tangent 
plane Tj.(S^). 

3. Lower bound for EA^(/) 

Here, we show that E,N{f) > n^. The proof has two ingredients: 
an estimate of the maximum max |/| and existence of the "barrier 

function" b^. 

Claim 3.1 (Estimate of the maximum). Given p > 0, there exists Cq 
such that, for any x G P{ max |/| ^ Cq} ^ |. 

Proof of Glaim \3.1[ By the mean- value inequality in Claim 12.21 for any 
spherical harmonic / G Tin and any x G we have 

max f^<n^ [ f. 

V(x,p/n) JV(x,{p+l)/n) 

Integrating this inequality with respect to x over the sphere, changing 
the integration order, and taking into account that n?'a{V{ ■ , ^^)) ^ 1, 
we get 

max /'< /" 
Hence, 

E{ max /^|=e| [ max f^] <E\\ff = 1 

(in the first equation we used the rotation invariance of the distribution 
of /). Applying Chebyshev's inequality, we get the estimate. □ 
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Claim 3.2 (Existence of tlie barrier). There exist positive numerical 
constants p and Ci such that, for each sufficiently large n and each 
X G S^, there is a function G Tin with the following properties: 

\%\\ = 1, b^{x) ^ CiVn, and ^ -CiVn. 

Proof of Claim ISTB: If x is the North Pole, then by Claims [231 12.81 and 
12.91 the zonal spherical harmonic Yq gives us what we need. For other 
x's, we just rotate the sphere. □ 



Proof of the lower bound for N{f) : Fix a; G S^. We have / = S^ob^ + fx 
where is a Gaussian random variable with E^q = 2;^, and fx is 
a Gaussian spherical harmonic built over the orthogonal complement 
to bx in Tin and normalized by E||/a.|p = i^^- We choose a Gaussian 

random variable ^0 independent of ^0 and of fx with E^q = E^q = 2;^, 

and set f± = ±^obx + fx- These are Gaussian spherical harmonics 
having the same distribution as /. Note that 

f = ^obx + \{f+ + f-) , 

and that by Claim 13.11 

P({max ^ Co} n {max |/_| ^ Co}) > 1 - {]- + h = \ . 

Now, consider the event Vtx that f{x) ^ Cq and /Igj,^^,.^/^-) ^ —Cq. 
The event VLx happens provided that 

io^/n ^ 2ci^Cq and max |/-|-| ^ Co . 

Therefore, 

P(n,) ^ P(eov^ ^ 2cr^Co) ■ 

■P({max ^ Co} n {max |/_| ^ Co}) ^ > . 

V{p,x) 'D(p,x) 

Here, k is a positive numerical constant. (Recall that the variance of 
the Gaussian random variable ^Qy/n is of constant size.) 

It remains to choose ~ disjoint disks on §^ of radius 2p/n. Each of 
them contains a component of Z{f) with probability at least n. Hence, 
EiV(/) > n\ □ 
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4. Exponential concentration near the median 

4.1. Main lemma. We would like to use Levy's concentration of mea- 
sure principle. To this end, we need to show that the number A^(/) 
doesn't change too much under slight perturbations of /. We won't 
be able to prove it for all / G TCn but we will show that the "unsta- 
ble" spherical harmonics / for which small perturbations can lead to 
a drastic decrease in the number of nodal lines are exponentially rare. 
More precisely, we will prove the following 

Lemma 4.1. For every e > 0, there exists p > and an exceptional set 
E cHn ofprobabiUty F{E) ^ C{e)e-''^^^'' such that for all f eUnXE 
and for all g G Tin satisfying \\g\\ ^ p, we have N{f + g) ^ N{f) —en'^ . 

Let us show that Lemma HIT] ensures an exponential concentration of 

— ^ near its median a„. Consider first the set F = {/ G 7i„: A^(/) > 

{an + e)n'^}. Then for / G (F \ we have A^(/) > a„n^, and 

therefore, P((F \ E)+p) ^ |. Hence, P(F \ E) ^ 26"'=^"" and 

Now consider the set G = {/ G 7Y„: N{f) < {an - e)n'^}. Then 
G+p cifeHn-. N{f) <al}UE 

and, thereby, 

P(G+,) ^ i + C(e)e-(^)" < ^ 

for large n and it follows that P(G) ^ 2e~'^''^" for large n. It remains 
to note that, for fixed c{e), we can always make the estimate hold for 
small n by increasing the value of C{6). □ 

4.2. Unstable spherical harmonics are exponentially rare. The 

exceptional set E of "unstable spherical harmonics" is constructed as 
follows. We take a sufficiently large positive R and cover the unit 
sphere by approximately R~^n? spherical disks T>j of radii R/n with 
multiplicity of covering bounded by a positive numerical constant. Let 
2)Dj be the disks of radii 3R/n with the same centers as T>j. Fix some 
small a, (3 > 0. We shall call a disk 3Vj stable for a function / G Tin 
if there is no point x G 3Vj such that |/(a;)| < a and |V/(a;)| < l3n 
simultaneously. Otherwise we shall call the disk 3Vj unstable. Finally, 
fix a small 6 > 0. We shall call a function / G 7i„ exceptional if the 
number of the unstable disks for this function exceeds Sn"^. 

Our first task will be to find the conditions that would imply that 
the exceptional functions are exponentially rare. To this end, note that 
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if we can find dn^ unstable disks, we can also find C25n^ unstable disks 
that are 4/n-separated. Now, for each unstable disk 3Pj in this well- 
separated family, pick a point Xj G 'iVj where \f\<a and |V/| < jSn 
simultaneously. Fix 7 G (0, 1) and consider the disks V{xj,'-y/n). They 
are pairwise disjoint. Let Mj = max |VV/|. Note that 

[ f> n-'Mf 

Jv{xj,2/n) 

and that the disks V{xj,2/n) are also pairwise disjoint. Hence, 

3 

Now there are 2 possibilities: either ||/|| > 2, or for the majority of our 
disks ZVj, we have ^ Ca^^^/^n^. 

The functions for which the first possibility holds are exponentially 
rare (Claim EHO]). 

On the other hand, if the second possibility holds, we can conclude 
using the Taylor formula that in at least |c25n^ pairwise disjoint disks 
of radius 7/^, we have the estimates 

I/I ^a + /?7 + C^2ri/V 
|V/| ^ {p^C^b~'l^^)n 

Let now g G 7Y„ satisfy \g\ ^ r. Then the number of our disks where 
max \g\ is much greater than b^^l'^T or max | Vf^l is much greater than 
b~^l'^Tn is small compared to br? . Thus, we can conclude that /+(? G U 
where U is the set of all h G 7Y„ satisfying 

Aih) = Area{x G §^ : \h{x)\ ^ A, \\/h{x)\ ^ Bn} ^ c^b-f^ . 

with 

A = a + + Cs6-'/\-f^ + t) 

B = P + C,6-'/\j + t) 

We want to show that F{U) ^ | and use Levy's concentration of mea- 
sure principle to conclude that the probability that / is exceptional does 
not exceed 2e~'^'^ By independence of h{x) and Vh{x) (Claim [2.12p . 
we see that, for each x G S^, we have 

¥{\h{x)\ ^ A, \Vh{x)\ ^ Bn} ^ dAB^ . 

Due to rotation invariance, 

EA{h) ^ F{\h{x)\ ^ A, \Vh{x)\ ^ Bn} ^ CiAB^ , 
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and we can draw the desired conclusion if 2C4,AB^ ^ 0357^. At this 
point we shall just note that, for given 5 > 0, we can always choose 
some positive 7, a, /3 and r to satisfy this inequality just because the 
right hand side behaves like 7^ and the left hand side behaves like 7^ 
when a = j3 = T = Q and 7 — 0+. We shall postpone the optimal 
choice of parameters until later when all the relations between them 
will be discerned. □ 

4.3. Now our task is to find the conditions that will ensure that N{f + 
g) ^ N{f) — £71^ whenever / is not exceptional and \\g\\ ^ p. We need 
to estimate the number of components of Z{f) that may disappear or 
merge with some other components in the process of perturbing / by 
9- 

First of all, we discard all components of Z{f) whose diameters is 
greater than R/n. Since the total length of Z{f) does not exceed Cn 
(Claim [2T3I) . we can conclude that the number of such components is 
much less than eri^ if R is much greater than e~^. 

Now, for each small component F, we fix the disk Vj that intersects 
F. Then F lies deeply within the disk "iVj-. the distance from F to the 
boundary of SPj is at least R/n. 

Next, we forget about all small components whose disks are unstable. 
The area estimate (Claim [231) implies that each unstable disk SPj can 
contain at most CR^ small components, so, if / is not exceptional, the 
total number of small components whose disks are unstable does not 
exceed 5i?^n^, which is much smaller than en^ if 5R^ is much less than 
e. 

We need to show that if the disk corresponding to the component F 
is stable, then the component F won't disappear or merge with another 
component unless max|5f| ^ a. This will follow from the next claim 

which we will use later in various contexts. 

Claim 4.2. Fix positive fi and v. Let D he a disk and let F he a 
-function on 2) such that at each point x G 2) either \F{x)\ > fi 
or |V/(a:)| > u. Then each component F of the zero set Z{F) with 
dist(T,dD) > ji/u is contained in an "annulus" Ar hounded hy two 
smooth curves and such that 

(i) F = +yU on one houndary curve of Ar and = — /i on the other; 
(li) Ar C F+^/^; 

(Hi) the annuli Ar are pairwise disjoint. 

Note that if G is an arbitrary continuous function on D with sup IG] < 
fi, then At must contain at least one component of the zero set Z{F + 
G). We get 
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Corollary 4.3. In the assumptions of the previous claim, suppose that 
G G C{D) with sup IGI < /i. Then each component T of Z{F) with 
dist (r, dD) > fJ'/iy generates a component T of the zero set Z{F + G) 
such that r C r^^/,^. Different components Ti ^ r2 of Z{F) generate 
different components Ti ^ of Z{F + G). 

Later we'll use this corollary in various contexts. 

Proof of Claim Replacing the function F{u) by AiF(A2u), we may 
assume that = z/ = 1. This will simplify our notation. 

Let us look at what happens with the connected component J^{t) of 
the set {\F\ < t} containing F as t increases from to 1. As long as 
jF(t) stays away from the boundary dD, it cannot merge with another 
component of {\F\ < t} because such a merge can occur only at a 
critical point of F and all critical values of F in X) are greater than t in 
absolute value. For the same reason neither of the two boundary curves 
of can collapse and disappear. But cannot come too close 
to dD before it merges with some other component either: indeed, if 
X e T{t) and T{t) hes at a positive distance from the boundary dD 
then we can go from x in the direction of V-F if F{x) < and in the 
direction —VF if F{x) > 0. In any case, since |V-F| > 1 in J-'{t), 
we shall reach the zero set Z{F) after going the unit length or less. 
Since the only component of Z{F) in J-'{t) before any merges is F, we 
conclude that J^(t) C F+i. Recalling that dist (F, dD) > 1, we see that, 
for each t ^ 1, J-'{t) stays away from the boundary dD. 

Thus, each component F lies in an "annulus" = J-'{1) which is 
contained with its boundary in the open disk D and such that F = 1 in 
one boundary curve of and F = — 1 on the other. By construction, 
the annuli Ar are pairwise disjoint. This proves the claim. □ 

Now, we apply Corollary 14.31 to the functions F = f and G = g on 
the disk D = 3Vj with fi = a, and u = j3n. We require that a/ 13 < R. 
This guarantees that if F is a component of Z{f) with diam (F) ^ R/n 
and T nVj ^ 0, then dist {T,d{3Vj)) ^ R/n > a/{(3n). We see 
that the only small components of Z{f) in stable disks Vj that can 
be destroyed by perturbation oi f hy g are those that correspond to 
the disks where max \g\ ^ a. By the mean value property (Claim [?^ . 

the number of such disks does not exceed Gp^a'^n"^ and, by the area 
estimate (Claim [23]), the number of the corresponding components is 
bounded by G p'^a~'^R^n^ , which is much less than en^ if is much 
less than ea^R"'^. 
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4.4. Tuning the parameters. Now it is time to malce tlie clioice of 
our parameters. First, let us list the constraints introduced above: 

< ea'^R'^, a < R(3, 5R^ < £, > 

and 

{a + (3^ + r + r)) {(3 + r ^^(^ + r))' « 5^^ . 
We take 

R ~ e~^, S ~ ~ a^e^, and /3 ~ ea. 

The quantity we want to maximize is c{£) ~ min(p^, r^) ~ min(T^, a^s^) 
subject to the constraint 

[a + e-'/'i^' + r)] ■ [e'a' + £-=^(7' + r')] « s'^' 

(we neglected absolute constants and the term /?7 ~ ae^y < a in the 
first bracket). Denoting the minimum to maximize by m, we see that 
we have to put r = m^/^, a — m^/^£~^/^. This leads to the constraint 

£-3/2[^l/2 ^ ^2] . ^-3[^2 ^ ^] ^< ^3^2 _ 

Again, we neglected e'^a'^ — e~^m < e~^m. Rewrite this constraint as 

[^1/2 + ^2] . ^ ^2] ^ ^15/2^2 _ 

It is immediate from here that m <S s^^ . On the other hand, taking 
^2 ^ £i5/2^ ^j^g^^ |-]^jg upper bound can be attained. Thus, the 

proof we presented gives c{e) ~ £^^. □ 

5. Existence of the limit lim a„ 

n— »oo 

In this section, we denote the spherical harmonics from Tin by /„. 
Since the random variable N[fn)/n^ exponentially concentrates near 
its median a„ and is uniformly bounded, it suffices to show that the 
sequence of means {EA^(/„)/n^} converges. Then the sequence of me- 
dians {a„} converges to the same limit. In what follows, we'll show 
that {EA'"(/„)/n^} is a Cauchy's sequence. 

5.1. Some integral geometry. Let ^ be a system of N{Q) loops on 
the sphere S^. By N^{Q^ T>) we denote the number of loops from Q that 
are contained in the spherical disk V, and by N{Q, V) we denote the 
number of loops from Q that intersect T>. We fix p > and denote 
T^x — T^{x, p), S — a{T>x). Note that the area S does not depend on x. 

Claim 5.1. 

If N4g,V,)da{x)^N{g)^l [ N{g,V,)da{x). 
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Proof: Fix a loop T G G and note that 

a{{x:T CV^}) S and S ^ a{{x: T nV,^ ^ 0}) . 

To prove the first inequahty, we fix an arbitrary point ?/ G F and observe 
that {x: F C V^} C Vy. Similarly, to prove the second inequality 
holds, we fix a point y G F and note that {x: T nV^ ^ 0} D Vy. □ 

Now, we fix 1 ^ (i ^ -R, put p = R/n, and let n go to oo. We set 

N{U,V,)''^ N{Z{U),V,) and N,.{U,V,) N,{ZUn),V,) . 

We call the component F of Z{fn) d-normalii its diameter does not ex- 
ceed c?, and denote by A^d(/n, ^^x) the number of (i-normal components 
of Z{fn) that intersect the disk T>x. 

By Claim 12.31 the total number of d-abnormal components does not 
exceed Crt^/d. Thus, applying Claim [5?T] to a spherical harmonic /„, 
we get 



1 



§2 



N,{fn,V^)d(j{x) 



71^ d 



^ j^^Ndifn,V,)da{x) + ^. 



Taking the expectation and using rotation invariance of the distribu- 
tion of random spherical harmonics (and recalling that cr(S^) = 1), we 
continue our chain of estimates 

EiV4/„,P,J ^ EiV(/„) ^ EiV,(/„,P,J ^ C 
Sn? ^ ^ Sn?' d 

where xq is an arbitrary point on S 



2 



5.2. Scaling. We fix a point xq G S^, denote by Xg the antipodal point, 
and by 7r^.„ : \ {xq} T^qS^ the stereographic projection (7r^o(xo) = 

0), and define a function F„ on T^,,S^ by F„,(u) = {fn°T^xl){~) ■ We also 

set D{t) = -D(0, t). By Nd{Fn, R) we denote the number of components 
of the nodal set {F„ = 0} of diameter at most d that intersect the disk 
D{R), and by N^{Fn,R) we denote the number of components that 
are contained in the disk D{R). Note that nxoi'Dx^) = D{0,R/n) with 
some R> R. Then iV^(F„,^) = A^^(/„, Since 1 < c? < < n, 
we have iVd(/„,P,J ^ N^diF^.R), 

^ R ^ ^(RV 1 . ^(RV Sn^ ^ 
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Then, using that KNd^fn, T^xo) ^ S, we get a scaled version of (15.21) : 
EiV,(F^,^) ^ ENjU) ^ EN2d{K,R) ^ CR^ ^ C 

vahd for 1 <^ d <^ R <^ n. At this point we simphfy our notation 
returning to notation d instead of 2d and to R instead of R. We get 

Claim 5.3. For any d and R such that 1 d <^ R <^ n, we have 

EiV(/,) EiV(/J ^ E{N,{F^,R)-N,{F^,R)} ^ CR^ ^ C 
V? 711? ^ i?^ r? d 

Later, estimating the expectation on the right-hand side, we'll use 

that the e.pte=.ion ^^l^- ^' " i= bounded from above by 

R^ 

a positive numerical constant. This follows from 



Claim 5.4. We have 

Nd{Fn,R) 
uniformly with respect to R and n. 



Proof: Obviously, Nd{Fn, R) ^ N^{Fn,R + d). By scaling Claim | 
the area of each nodal domain of the function F„ cannot be less than a 
positive numerical constant C4, therefore, A^*(F„, R+d) ^ 7r(i?+(i)^cJ^. 
Hence, the claim. □ 

In what follows, we show that if we discard some events of small 
probability, the difference Nd{Fn,R) — N^:{Fm,R) will be small. In 
view of Claims [5?3] and [5^ this will prove that EN{fn) /n^ is a Cauchy's 
sequence. 

The main idea is to show first that if m and n are sufficiently 
large, then the function Fm can be viewed as a statistically small C^- 
perturbation of the function F^, and therefore, outside of small events, 
Nd{Fn, R) cannot be much larger than N^{Fm, R). We start with 

Claim 5.5. Given a finite set of points {uj^ , the random vectors 
Fn{uj) converge in distribution as n ^ 00. 

Proof: We use Claim O 

1 " 

E{/n(x)/„(2/)} = M^)^k{y) = Pn{cOSe{x,y)) 

k=—n 



ON THE NUMBER OF NODAL DOMAINS 15 

where 9(x, y) is the angle between x and y as vectors in M^. Then the 
scaled covariance equals 

E{F4u)F„(v)} = e((/„ oO(^)(/„ o^-;)(^) 

When n goes to oo, the angle between the points 7i"xoUn)' ''^xoin) 
on the sphere is equivalent to \u — v\/n (locally uniformly in u and 
v). Therefore, by Hilb's theorem (Claim I2.9p . the scaled covariance 
converges to the Bessel kernel Jod'^ ~ ^1) locally uni- 
formly in u and v. 

Recall that the vector |F„(-Uj)} is a Gaussian one, and that the 
convergence of covariance matrices of a sequence of Gaussian vectors 
yields convergence in distribution of the vectors. □ 

5.3. Discarding small events. 

5.3.1. Consider the event 

= F^dm>R']. 

Since at any point x G E|/„(x)p = 1, we have 

E I F^dm= I EF^ dm = CR^ . 

Jd{5R) Jd{5R) 

Then, by Chebyshev's inequality, F{n^n^) < R-^ and P(fii^^ U fi^^^) < 
R~^. Throwing away these events, we assume that 



2 / 172 
n 5 / 

D(5R) Jd{5R) 



max <! / F^, / F^[ ^R\ 



By Claim 12.21 this yields the estimates 

(5.6) ||-^n||c2{D(4i?)) , ||-^m||c2(D{4K)) ^ -R^^^ • 

5.3.2. Now, we fix a finite i?-('^+2)-net {uj} in the disk D{AR). The 
parameter a > 1 will be chosen later. Since by Claim 15.51 



lim P( max I (F„ - F^) {uj) | > e) = , 

min(m,n)^oo j 



in what follows, we discard the event 
nil 



I max I (F„ - FJ (u,) | > } 



16 FEDOR NAZAROV AND MIKHAIL SODIN 

and assume that 

max I (F„ - [uj) \ ^ . 
Using a priori estimates (15.61) . we get 

max \F„ — Fm\ ^ „ , „ H — „ , ^ < 



Then, scahng local gradient estimates from Claim [2^21 we get 

C 1 
max |V(F„ - ^— tt^ < ^ 

D{3R) ' ' 

if R is big enough. We conclude that 

1 

that is, outside of events fin'' U and fii^L, the random function 
Fjn indeed can be viewed as a small C^-perturbation of the random 
function F„ in the disk D{3R). 

5.3.3. To be sure that a i?~ "-perturbation of the function F„ does 
not decrease drastically the number of the components of the zero set 
|F„ = O} in the disk D{R), we need to know that the function F„ is 
"stable" in a larger disk D{3R), e.g., that 

min {|F„| + |VF„|} > • 

D{3R) ' ' 'J Ra 

Suppose that 

min {\Fn\ + |VF„|} ^ J- , 



and estimate the probability of this event (we call it fin )• We fix a 
^-(a+2)_j^g^ {mj}, this time in the disk D{3R), that contains at most 
CR^'^+e elements. 

Suppose that, at some point u G D{3R), 

|Fn(M)| + |VFn(w)|^ J^. 

Then there is a point Ui of our net such that 

(5.7) |F„(n,)| + |VFn(n,)| ^ + < ^ 

(we again used a priori estimates (15. 6p ). 

By the independence Claim 12.121 the probability that in a given 
point Uj from our net condition (15.71) holds does not exceed CR"^"^. 
Hence, the probability that (15.71) holds at some point of the net does 
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not exceed CR^'^^^-CR = CR^ ^ and tends to as i? — > oo provided 
that a > 6. Hence, choosing a = 7, we achieve that P(il^f^) ^ CR~^. 

To summarize, we denote by Q* the complement to the union of our 
small events n^n^ U fi^^^ U U fi^f^. Then 

C 

¥{Q* doesn't occur ) ^ h /t(n, m) , lim K,{n, m) = . 

min(n,m)— >oo 

We have proved the following 

Claim 5.8. Given R, n and m such that 1 <^ i? <^ min(n, m), there 
exists an event Q* such that if it happens then 

\\Fn — Fm\\c^(D(3R)) ^ , 



mm 



2 



i:)(3R) ^' ' ' 'J _R7 



and 

EN{fn) EN{f„ 



with 



, iVrf(F„,i?)-iV,(F„,i?) i?2 1 
< sup + — + 3 + '^l^' 

r2* occurs -fL TL d 



lim K(n, m) = 0, and 1 <^ d <^ R. 

min(n,m)— »oo 



5.4. The following claim estimates the supremum on the right-hand 
side of the previous bound. 

Claim 5.9. // the event Vt* occurs, then 

N,{F^,R)-N,{F^,R) ^ Cd 

Proof: First, note that Nd{Fn,R) — N^{Fjn,R) does not exceed the 
number of components of the zero set {Fm = 0} that are contained 
in the annulus A = |i? — c? — 2 ^ |m| ^ i? + + 2}. To see this, we 
apply Corollary 14.31 to the functions F = Fn and G = Fm — Fn with 
D = D{3R), and ^ = v = R~^ . By this Corollary, each ci-normal 
component F of Z{Fn) such that rr\D{R) ^ generates a component 
F of the zero set Z{Fm) such that F C F+i. If F is not contained in 
D{R) then, observing that the diameter of F does not exceed d + 2,we 
conclude that it must be contained in the annulus A. 

Since the area of each nodal domain of Fm cannot be less than a 
positive numerical constant (Claim 12.41) . we see that the number of 
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components of the zero set {F^ = 0} that are contained in the annulus 
A cannot exceed 

C[{R + d + 2f -{R-d-2f)] ^CdR, 

proving the claim. □ 

Combining Claims [518] and \5M we obtain 

2 2 ~ D + 3 + — + 0^' ^) ' 

First, we set d = \/R. Then, given £ > 0, we choose R so big that 
1/\/R < e. At last, we choose n and m so large that ^(n, m) < e and 
R'^/n^ < e. Then we get 

EiV(/„) EN{f„ 



71^ im? 



< 



This completes the proof of convergence of ]EN{fn)/n'^, and hence fin- 
ishes off the proof of the theorem. □ 

6. Sharpness of Theorem [TI} P{A^(/) < nn^} ^ e^^^^^)" 

The idea is very simple: the zero set of the zonal spherical harmonic 
Yq is a union of n circles of constant latitude. On the other hand, by 
Claim 12. 6[ the zonal harmonic Yq is stable, and therefore, its small L^- 
perturbations cannot increase much the number of components of the 
nodal set. 

Let 

n 

k=—n 

be a Gaussian spherical harmonic of degree n. Consider the event 

where p is a small positive constant which we shall choose later. We 
have 

In what follows, we assume that the event Q occurs. Then / = ^oYq + q 
with ll^fll ^ p. 

Again, we cover the sphere by ^ R~'^n^ spherical disks Vj of 
radius R/n with R (depending on n) to be chosen later. The disk 
Vj is good if irm,x (^\g\ + ^iVf?!) < \cq where cq such a constant that 
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(^0^0)^ + ^|V(^o>o)P ^ Cn everywhere on (ClaimEjl). Comparing 

the areas ( Claim [?]2l) . we see that the number of had disks is < c^'^p^n'^. 

Now, let r be a connected component of the nodal set Z{f). Then 
at least one of the three following possibilities must occur: 

(i) The component F has diameter larger than R/n. 

(ii) The component F has diameter less than R/n and intersects a good 
disk Vj. 

(iii) The component F has diameter less than R/n and intersects a bad 
disk Vj. 

By Claim [2731 the number of components of the first type is bounded 
by CR~^n^ -C /tn^, provided that R is chosen sufficiently large. 

If the disk Vj is good, then, for x G SPj, we cannot have < 
\cq and |V/(x)| < jCq at the same time. Therefore, we can apply 
Corollary 14.31 to the functions F = f and G = —g with = |co and 
u = |co. We see that each component F of the zero set Z{f) of the 

second type generates a component F of the zero set Z{Yq) of diameter 
at most {R + 2)/n. Recall that Yq is a zonal spherical harmonic and its 
nodal set consists of the spherical circumferences of constant latitude 
that are generated by zeroes of Legendre polynomials. By Claim 12.81 
the components F of diameter at most {R + 2)/n must be located 
in a neighbourhood of one of the Poles, and there is only a bounded 
number of them. Hence, the number of components F of the second 
type remains bounded as n goes to 00. 

At last, all components of the third type are contained in the set 
3Vj of area < p^R^, and by the area estimate (Claim [53]) the 

'Dj is bad 

number of such components is < p^R^n^ <^ nn"^ provided that p is 
properly chosen. □ 
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